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Abstract 

It is discussed the limitations of the widely used markovian approximation applied 
to model the turbulent refractive index in lightwave propagation. It is well-known 
the index is a passive scalar field. Thus, the actual knowledge about these quantities 
is used to propose an alternative stochastic process to the markovian approximation: 
the fractional Brownian motion. This generalizes the former introducing memory; 
that is, there is correlation along the propagation path. 
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1 Introduction 



Whenever a light beam propagates through the turbulent atmosphere, it ex- 
periments deflections due to fluctuations in the refractive index. As a result 
of phase changes, the beam suffers displacements perpendicular to the origi- 
nal direction of propagation. This phenomenon is known as beam wandering, 
beam steering or spot dancing. The wandering is usually characterized in terms 
of its variance. Several authors have experimentally and theoretically treated 
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this problem, or the equivalent problem of angular fluctuations, using different 
approaches. 

The earliest studies of the problem of beam wandering were within the Geo- 
metric Optics approximation. Chernov treated the ray-light propagation in 
a random medium as a continuous Markov process. This assumption enabled 
him to formulate a Fokker-Planck equation. Later, Beckmann j.2| used the ray 
equation to calculate the wander of a single ray. Both authors obtained sim- 
ilar formulas for the variance of the transverse displacements, Q = {Qx,Qy), 
observed after the light has propagated a distance L. These formulas differ in 
a numerical factor, but the power-law 

VarQocL^ (1) 

is found in both. 

However, Chernov and Beckmann introduced an arbitrary Gaussian covariance 
function for the refractive index fluctuation which is not physically plausible; 
also, these techniques do not include the effects of a flnite beam diameter. 
That is, small-scale variations of the refractive index contribute to the beam 
spread but have little effect on the motion of the beam centroid. Therefore, 
Geometric Optics approximation tends to overestimate the magnitude of the 
wandering. Chiba [sf tried to overcome this limitation. He assumed that only 
changes of the refractive index larger than beam diameter contribute to the 
wandering. Using a Kolmogorov-like structure function he also found a power- 
law dependence as in Eq. (1). 

On the other hand, whenever the width of the beam is smaller than the tur- 
bulence inner scale the approximations obtained from Geometric Optics are 
enough. These beams known as thin-beam were studied by Consortini and 
O'Donnell j^. Following the Beckmann paper, they analized experimentally 
and theoretically the dependence of thin-beam displacements with the propa- 
gation distance. They showed that, independently of the turbulence spectrum, 
the displacements variance grows like the third-power of the path length. That 
is, provided the fluctuations are small and the propagation length L is large 
compared against the outer scale Lq. In another paper, Consortini et al., |5| 
investigated experimentally, for laboratory-generated turbulence, the depen- 
dence of beam variance on propagation length in the case of strong turbulence. 
They also found a rough estimate of the third-power law dependence. 

It must be stressed here that all the previous works are valid in the Geometric 
Optics limit (see Ref. 0, p. 120). That is, 

lo > ^/LX, 

where A is the wavelength and Iq the dimension of the inner scale. 
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Klyatskin and Kon 0] described the propagation of light considering a scalar 
parabolic equation and assuming a markovian approximation for the atmo- 
spheric refractive index fluctuation, e(p; z) — p represents transversal coordi- 
nates (x, y). Using this approximation they derived an equation for the angular 
beam wandering, 6 = dQ/dz. They found a dependence proportional to L in 
accordance with the results above — the displacement and angular beam wan- 
dering variances are proportional with factor (see Ref. jj], Eq. (15)). 

When using the markovian approximation the original covariance function R, 
obtained from the original structure function for the index, is replaced by an 
effective covariance A through: 

E[e{p,z)e{p',z')]=6{z-z')A{p-p'), (2) 

where 



A{p) = I R{p,z)dz 

JR 



is a differentiable function, and E[-] is taken with a gaussian probability mea- 
sure. This means that the values of the index in the region z > z' do not 
affect those at the point (p', z') — which eliminates backscattering and scatter- 
ing at sufficiently large angles. This property known as "dynamical causality" 
(see Ref. j|8|, p. 214) comes from the martingale property of the markovian 
approximation. ^ 



Tatarski and Zavorotny j8| derived the conditions of validity for the marko- 
vian approximation. They found that it is applicable if all the characteristic 
dimensions arising from the wave propagation problem are small compared to 
the path length. As Ostoja-Starzewski explained jl3] an intuitive justification 
for the Markov property is that the ray-light on a long distance behaves as if 
it has suffered many independent refractions. Then, this approach holds only 
for long-path propagation. 

In Sec. 2 the equivalence between the markovian approximation and the pro- 
cedures usually followed after the use of the Geometric Optics approximation, 
referred above, is shown. Thus, the former model is only valid for long-path 
propagation and weak turbulence, or just strong turbulence. 

Recently, it was shown that fractional Brownian motion processes (fBm) could 



be used to describe the turbulent refractive index fluctuation [ll[ to model ray- 
light propagation. This is not new since these fractal stochastic processes have 
also been used to identify turbulence degraded wave- fronts [llEl- In Sec. 3 



we will show that these processes match the requirements for passive scalar 



^ Formally, it is said that a process Xf is markovian or possesses the Markov prop- 
erty if the future behavior of it given what has happened up to time t is the same 
as the behavior obtaiiied when starting the process at Xt — a detailed description 
can be found in Ref. 
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fields; that is, they are nearly Gaussian, have stationary increments, and obey 
the Kolmogorov-Obukhov-Corrsin (KOC) structure function — the extension to 
passive scalars of the well-known Kolmogorov structure function. Moreover, 



it was verified through Wavelet Analysis of experimental data, that the 
wandering of a laser beam presents memory — for short distances. This is in 
accordance with the fBm model presented above. Otherwise, the markovian 
approximation is memoryless; therefore, it can not represent the behavior of 
the refractive index in all spatial scales. Under these circumstances the fBm 
is a good candidate to replace and extend the markovian approximation. 



2 Markovian approximation 



Most of the past and present research in turbulent lightwave propagation, 
we have seen, is based directly or indirectly on the markovian approximation. 
Usually, when approximations from the Geometric Optics are used the marko- 
vian property is indirectly applied [H Q, 0, 0| • These papers propose a simple 
model: a ray-light beam (laser beam) propagates through a turbulent flow a 
distance L and the position on a screen is evaluated. 

The Geometric Optics ray-equation under the condition e <^ 1, small index 
pertubation, is 

= Vq6(Q,z). (3) 

The assumption leads to two independent differential equations; therefore, it 
is enough here to consider just one axis. Usually, what the literature refers as 
Geometric Optics approximation is obtained taking the zeroth order approx- 
imation for the transverse displacement, Qx, from the above equation. That 
is, 

rL rL rz rz' Q'^B^ 

VarQ^=/ / / / \ {p,s; p,s')dsds'dzdz', (4) 

JO Jo Jo Jo oxox 

where i?e(r, r') denotes the covariance function of the index fluctuation at 
points r = (p, s) and r' = {p',s'). Observe that to obtain this equation it 
should be assumed that the gradient of the turbulent index is bounded. 

It is always supposed the turbulence is homogeneous and isotropic. Under 
these hypotheses the covariance of the index fluctuation, B^, is assumed to be 
an even function of the coordinate difference: 

5,(r,r')=5.(||r-r'||), (5) 

where r = {p,z), r' = {p',z') and, as usual, || ■ || is the norm. After some 
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elaborated calculations the following result is obtained 







f{L,z)^iO,z)dz. 



(6) 



where f{L, z) = {2L^ — 3L^z + z^)/3 is called filter function. All the papers 
considered that the covariance has only significant values for z <^ L. Therefore, 
the second and third terms in the filter function can be neglected changing 
the latter equation to: 



VarQa 



2L3 fLd^B, 



(0, z) dz 



2L^ 



d'^x 



(0, z) dz. 



(7) 



This result holds for L ^ Lq, long-path propagation. 



Approximations involved in the previous relation are equivalent to the marko- 
vian approximation. This can be checked as follows. Let us introduce the 
covariance function (2) in Eq. (4): 



Var 



"'0 JO JO 



dx' 



dx"^ 

rL 



5{s — s)ds ds'dz dz' 



p=0 



p=0 



f ds f dzQ{z-s') I dz'e{z' 
Jo Jo JO 



dx"^ 



p=0 



3 ■ 



It is important to note that the appearance of a third-power law is character- 
istic from the markovian approximation. 

Since, it is considered an isotropic turbulence the effective covariance function 
is written as (see Ref. [l^, §20.2), 



A{p) = (27r) 



K,dKJo{Kp)^{K). 



(9) 



Thus, substituting the equation above in Eq. (8) it is exactly reproduced the 
result obtained by Consortini and O'Donnell for the long-path case: 



47r2 

VarQ,. = — 



up 



(10) 



p=0 



It should be stressed that Consortini and O'Donnell obtained this result from 
Eq. (4), so they are within the Geometric Optics approximation. Therefore, 
it is verified that this approach uses indirectly the markovian approximation. 
Moreover, it is concluded that the markovian approximation is valid for long- 
path and weak turbulence as it was mentioned at the end of the introduction. 
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On the other hand, if the propagation length L is small the complete vari- 
ance, Eq. (6), should be calculated. In Ref. 0] it is numerically evaluated. 
The variance is not proportional to the third-power law of the lenght L, but 
asymptotically approaches to it for long-path propagation. 

Consortini et al. [1] also found a cubic dependence in the strong turbulence 
case. Having into account that, in this case, the refractive index covariance 
along the direction of propagation has little effect on the characteristic fluctu- 
ation of the ray, the markovian approximation should also be valid for strong 
turbulence. 

Finally, a simple dimensional argument can show the intrinsic third-power 
law behavior of the markovian approximation. Observe from Eq. (2) that the 
refractive index fluctuation behaves as a white noise W^^'^ along the z-axis, so 
Vqc oc W^I"^ . Henceforth, using Eq. (3), the angular beam wandering Q grows 
as L^/^. The variance of the displacements is proportional to . 



3 Fractional Brownian motion and the turbulent refractive index 

As it is widely known the atmospheric refractive index inherits from the tem- 
perature field the quality of being a passive scalar field [3l — whenever the 
temperature has also been proved to be a scalar field. Over the last decade the 
quest of the Fluid Dynamics has been to answer the following question: what 
properties are passed down by the stochastic nature of the turbulence to the 
scalars embedded within? The answer is not simple neither has been closed. 

Any stochastic variable is characterized through its probability distribution, 
so the scalars are. First note that it is always possible to assume that any 
scalar "d is locally homogeneous, just as the velocity field, and therefore its 
increments are (spatially) stationary random processes — see Ref. 0, p. 19. 
Thus, it implies that for any statistical moment 

(Kr)-^(r')n = ([t9(r-r')- W), 

where ( ■ ) stands for the average over an undefined probability. Nevertheless, it 
must be stressed that the scalar itself is not stationary; moreover, the gradient 
of the averaged scalar determines much of its statistical behavior. 

Therefore, it is preferred to study the moments of the increments: the n-point 
structure functions. It has been confirmed that these structure functions vanish 



^ Remember, that such quantities are diffused and advected by the turbulence with 
negUgible back effect on the flow. 
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if n is odd, and behave as follows otherwise 

52n(r) = ([^(r)-^(0)]2") = A„||r||f-, 



'111 



for ||r|| < Lq. The coefficients C2n are controlled by the velocity field, an 
(isotropic) external force acting over the scalar and the gradient of the mean 
scalar concentration. Given the coefficient (2, a Gaussian behaviour should 
be observed only when (2n = JK2- Unfortunately, there is plenty of evidence 



that this is not the case, e. g. [17[: these exponents deviate from the Gaussian 
estimated value. Also, the very same coefficient (2 depends on the intermit- 
tent behavior of the velocity field itself, and it is different from the '2/3' 
Kolmogorov exponent. It has been observed to range from 0, as the visc osity 
becomes relevant, to 2, as the external force reaches the innermost scales [l8|. 
An additive factor called intermittency exponent can be introduced to measure 
these deviations — see for example Ref. HI HIS. 



Nevertheless, only the first moments of the propagated light are relevant to 
us. For that, the Gaussian distribution should be enough. This is how it has 
been done while studying scalar turbulence from synthetic Gaussian velocity 
fields — see again references [l^ 12) 



Finally, the 2-point structure function, Eq. (11) with n = 1, is called KOC 
structure function or the so called structure function. It is related to the 
covariance function f(r, r') = {■}}{r)'d{r')): 



S2{r - r') = ^;(r, r) + w(r', r') - 2w(r, r'). 



(12) 



Observe that the covariance function is symmetric under the change r ^ r'; 
thus, dv{r,r)/dx = dv{r,r)/dx' — the same happens in the other axis. From 
Eq. (12) it is 



dS2 
dx 



dv ^ 
dx 



dv , 



Then, from the latter and Eq. (11), it is obtained: 



dv dv 
T^(r,r) - T^(r,r') 
dx dx 



C2A2 



[X — X 



r'l|2-C2 



(13) 



where 2 — (^2 > 0; therefore, as r ^ r' the left-hand side diverges. As it is 
known a process whose covariance function lacks second partial derivatives is 
not derivable (see Ref. i23|, §9.4); thus, a stochastic process meant to model 
a scalar should obey this property. 



On the other hand, Stolovitzky and Sreenivasan i^j successfully obtained the 
Kolmogorov's law — the velocity version of Eq. (11) — using fractional Brow- 
nian motion (fBm) to model the turbulent velocity field. Therefore, since 
the phenomenological parallelism between the turbulent velocity field and the 
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scalar fluctuations [13] these stochastic processes can also be considered here. 
It can be seen that they fulfill indeed all the properties mentioned before for 
the refractive index. 

Briefiy, fractional Brownian motions are a family of Gaussian processes , 
being H the Hurst parameter [25'], with covariance j26| : 



E 



= i (|t|'^ + \s\'" -\t- , (14) 



2 



for s,t G M, < if < 1, 5^(0) = almost surely, and E[5^(t)] = 0. Note 
that these processes are non-differentiable. Also, they are scalar- invariant; 
that is, B^{as) = a^B^{s), for any a, where = means both share the same 
probability law. Usually, scalar-invariant processes are called self-similar if 
they have stationary increments. This is effectively what happens here, as it 
can be observed from Eq. (14). 

One remarkable property of this family -B^ is that the H parameter regulates 
the presence or absence of memory [23|. In fact, it can be separated in three 
subfamilies accordingly: long- memory for 1/2 < if < 1, no-memory at ii = 
1/2, and short-memory in the case < ii < 1/2. 

Now, let the isotropic fractional Brownian motion (zfBm) be defined as 

B^{r) ■.= B^{\\r\\) = B"{r). 

It is straightforward to calculate the covariance of its increments from Eq. (14). 
Afterwards, it is observed that only when \r — r'f^'^ <^ 1, 



E 



[B"{r)-B 



~iir-r'f^. (15) 



Consequently, the turbulent refractive index can be modeled as follows 

e(r) :=a5^(r/Lo), (16) 

where a is an adimensional constant. Using the approximation in Eq. (15) it 
is straightforwardly obtained the structure function as below: 



E 



(e(r + r')-e(r))' 



a L, 



-2H\\„\\2H 



for ||r|| ^ Lq. Thus, comparing the latter against Eq. (11) for n = 1, it 
results C2 = 2ii (0 < ii < 1) and A2 = a^L^'^^ — using typical values of the 
structure constant for the refractive index and the outer scale a ~ 10^^-10^^. 

Again, a dimensional analysis on the Geometric Optics ray-equation, Eq. (3), 
can be made to find the path dependence of the refractive index. From Eq. (16), 
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for IIQII < 1, it is: 



WzBz + Qll Z 



where is the fBm noise. The noise, thought not a derivative in the usual 
sense, extends the chain rule so the above is valid. Since, Q/z 9 and a is 
small then 9 oc . Therefore, its variance should grow as L?^ and for the 
displacements as L?^^"^. 



4 Conclusions 



The procedures employed after the Geometric Optics approximation, which is 
usually used in atmospheric optics, has been shown equivalent to the marko- 
vian approximation. Here is verified its validity in the following cases: at prop- 
agation path long enough for weak turbulence, or at any path length for strong 
tubulence. Moreover, as the cubic dependence on the path length is inherent 
to the markovian approximation, the Geometric Optics approximations, Eqs. 
(5)-(7), are indirectly based on the former. But, it means that no-memory 
processes are related with. 

A new model for the turbulent refractive index fluctuations is proposed in this 
paper following most of the requirements of a passive scalar. The isotropic frac- 
tional Brownian motion extends the markovian approximation, and introduces 
memory in the ray-light propagation phenomenon. Furthermore, through the 
Hurst parameter the state of the scalar turbulence can be set. Then, the vari- 
ance of the displacements behaves as L^^+^. In particular, for Markov pro- 
cesses H = 1/2, and the cubic dependence is recovered. This is the situation 
of long path propagation for weak turbulence, and strong turbulence at 
any path length |5| discussed in Sec. 2. 

Finally, the presence of memory, H ^ 1/2, was experimentally observed for 



short path propagation jlj|. Nevertheless, the power-law dependence men- 
tioned above must be experimentally confirmed. That task will be the chal- 
lenge of future works. 
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